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Probing the LMC age gap at intermediate cluster masses 
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ABSTRACT 

The LMC has a rich star cluster system spanning a wide range of ages and masses. 
One striking feature of the LMC cluster system is the existence of an age gap between 
3-10 Gyrs. But this feature is not as clearly seen among field stars. Three LMC fields 
containing relatively poor and sparse clusters whose integrated colours are consistent 
with those of intermediate age simple stellar populations have been imaged in BVI 
with the Optical Imager (SOI) at the Southern Telescope for Astrophysical Research 
(SOAR). A total of 6 clusters, 5 of them with estimated initial masses M < 10* Mq, 
were studied in these fields. Photometry was performed and Colour-Magnitude Dia- 
grams (CMD) were built using standard point spread function fitting methods. The 
faintest stars measured reach V ^ 23. The CMD was cleaned from field contamination 
by making use of the three-dimensional colour and magnitude space available in order 
to select stars in excess relative to the field. A statistical CMD comparison method 
was developed for this purpose. The subtraction method has proven to be successful, 
yielding cleaned CMDs consistent with a simple stellar population. The intermedi- 
ate age candidates were found to be the oldest in our sample, with ages between 1-2 
Gyrs. The remaining clusters found in the SOAR/SOI have ages ranging from 100 to 
200 Myrs. Our analysis has conclusively shown that none of the relatively low-mass 
clusters studied by us belongs to the LMC age-gap. 

Key words: galaxies: Magellanic Clouds; galaxies: star clusters; galaxies: stellar 
content; stars: statistics 



INTRODUCTION 

' The Large and Small Magellanic Clouds (LMC and SMC, 
respectively) make up a very nearby system of low-mass, 
gas-rich and interacting gala^xies. At their distances, both 
Clouds can be resolved into stars, allowing their stellar pop- 
ulations and star formation history (SFH) to be studied 
in detail. These studies open up the possibility to identify 
epochs of enhanced or reduced star f ormation, and to asso- 
ciate them to the system dynam ics l|Holtzman et al.lll999l : 
iJaviel. Santiago, fc KerbedbOOSl ). 

Star clusters provide an alternative tool to reconstruct 
the SFH of a galaxy. The Magellanic Clouds have a large 
cluster system, spanning a wide range of properties, such as 
masses, ages and metallicities, which, gi ven their proxi mity, 
can be determined using different tools ('S antiagoll2009l ). Al- 
though they make up only a small fraction of the stellar mass 
in a given galaxy, star clusters have an advantage over field 
stars in that they may be modelled as simple stellar popula- 
tions (SSPs), facilitating derivation of their main properties. 
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Sizeable samples of clusters with available integrated mag- 
nitudes and colours currently exist. In the LMC, they have 
been used in association with evolutionary SSPs models to 
derive age and mass distributions and to reconstruct the 
Initial Cluster Mass Function (ICMF) and the Cluster For- 
mation Rate (CFR) jHunter et al.ll2003l : Ide Griis fc Anders! 
l2006l : IParmentier fc de Griisll2008l ). 

One striking and undisputable feature of the LMC clus- 
ter system is the so-c a lled a ge gap. It was proposed by 
Ijensen. Mould, fc ReidI (|l988l ). as a lack of clusters in the 
range 4 ^ t ^ 10 Gyrs. Several candidat es to fill this gap 
were proposed and discarded since then (|Saraiedinilll998l : 
Rich, Shara, fc Zurek 200lh. Evidence for the gap is presen t 
in the CFR reconstruction bv lParmentier fc de Grii3 l|2008f ). 
An age-gap is also tentatively seen in some reconstructions 
of the LMC SFH based on field stars, although not system- 
atically in all fi elds studied, and not as clearl y as in the 
clust er system (jjaviel. Santiago, fc Kerb er 2005.: [Noel et al.l 
l2007h . The lack of a clear age gap among field stars suggests 
that it may be less pronounced among lower mass clusters 
(M < IO^Mq) as well, which tend to be systematically un- 
favoured in current magnitude- limited cluster samples in the 
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LMC. Thus, the sample of candidate clusters with well de- 
termined ages must be pushed towards fainter limits than 
has been previously done. 

The most reliable cluster age determinations require 
CMD analysis. Confrontation of observed CMDs with the- 
oretical ones, built from stellar evolutionary models, pro- 
vide the best tool for age determinatio n of single clusters , 
either with HST or fro m the ground (|Baume et al.l l2008l : 
IPiatti et al.ll2009l . 12007( 1. Building a cluster CMD, however, 
is a more costly task than obtaining integrated photometry, 
as it requires larger telescopes with high pixel sampling and 
good seeing conditions. Therefore, selection of candidates to 
fill the LMC age gap should be based on the afore mentioned 
samples with integrated photometry. Furthermore, poor and 
sparse clusters suffer from more severe field star contamina- 
tion on their CMDs, something that leads to a bias towards 
rich LMC clusters having detailed CMD analysis available. 

In this paper we analyse the CMDs of 6 LMC clusters 
for which no CMD is yet available. Two of these clusters have 
integrated colours consistent w ith an intermediate age, ac- 
cording to the photometry from llfunter et al] (|2003l ). Five of 
them are sparse and unconspicuous when compared to pre- 
vious samples. Our main goal is to probe relatively poor and 
intermediate-age clusters in search of possible examples that 
may have ages within the age gap. They have been imaged 
with the Optical Imager (SOI) at the Southern Telescope 
for Astrophysical Research (SOAR), under 1" seeing and 
with excellent PSF sampling. In Sect. 2 we describe the ob- 
servations and photometry. In Sect. 3 we present our anal- 
ysis tools, including a method to efficiently decontaminate 
the cluster CMD from field stars. This new decontamination 
method uses the entire information on the magnitude and 
colour space to subtract field stars from the CMD. We show 
the results for each cluster individually. In Sect. 4 we present 
our final discussion and conclusions. 



2 DATA 

The observations were made in service mode and took place 
in the nights of November 10th and December 16th 2007, 
using the SOAR Optical Imager (SOI) in B, V and I filters. 
Fie lds were imag e d arou nd 3 LMC age-gap candidates listed 
by iHunter et al.1 (|2003D Q : OGLE-LMC0531, KMK88-38, 
and OGLE-LMC0169. Another field was observed around 
the richer and yet unstudied cluster NGC 1878. The field 
around OGLE-LMC0169 was observed under unstable and 
mostly bad seeing for our purposes, which led us to dis- 
card the data. For the remaining fields, the mean seeing was 
around 0.9 arcsec during both nights. 

The SOAR Optical Imager (SOI) uses two 2050 x 4100 
px CCDs, covering a 5.26 arcmin square field of view (FOV) 
at a scale of 0.077" /pixel. The images have a gap of 10.8" 
between the two CCDs. A slow readout was adopted in order 
to minimize detector noise. A 2x2 binning was used, yielding 
a spatial scale of 0.154" /pixel. 

In Table [1] we summarize the observation log, where 
we show: target name, filter, seeing, number of exposures. 
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Figure 1. Sections of the final V band images wliere the clusters 
arc located. The size of each section is 50" X 50". The cluster 
in each panel is as follows: (a) KMK88-38 (b) KMK88-39 (c) 
OGLE-LMC0531 (d) OGLE-LMC0523 (e) NGC 1878 (f) OGLE- 
LMC0214 .In all figures North is up and East is to the right. 



and exposure time. Multiple exposures were taken in or- 
der to increase the signal to noise ratio and reject cosmic 
rays in the final combined image. For each cluster, a set of 
short exposures was also taken in order to avoid saturation 
of bright stars. The seeing is the mean full width at half 
maximum (FWHM) of bright stars from the individual ex- 
posures. These exposures were flatfielded, bias subtracted, 
mosaiced and combined. 

The SOAR/SOI fields i ncluded 3 add itional star clus- 
ters besides those taken from lHunter et ahl (12003, ') . They are 
OGLE-LMC0214, OGLE-LMC 0523, and KMK88- 39. The 
clusters names are taken from iHunter et all (|2003l ). Fig. \T\ 
shows the sections (50"x50") of the final combined V band 
images where the clusters are located. Five of them are low 
density and relatively sparse open clusters; only NGC 1878 
is a richer and more compact object. To our knowledge, none 
of them have published CMD data in the literature. 

For photometric calibration, two standard fields from 
ISharpee et all (|2002l ). located at the north-eastern arm of 
the SMC, were observed on each night. The observations 
have distinct air masses ranging from 1.34 to 1.76. The mag- 
nitudes obtained from these images were used to fit a cali- 
bration equation for each passband (B,V, and I). We took 
care to include standards covering a similar airmass range 
as the clusters, and under seeing ^ l". The fit residuals were 
~ 0.03 mags, indicating that both nights were photometric 
for most of the time. 
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Table 1. Log of observations. 
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generating a master list of stars, (ii) This master list was 
then used as the position reference, (iii) The photometry 
on the remaining passbands then uses this master list, with 
only small positional offsets applied. 

In Fig. [2] we show the residuals between the photometry 
of long and short exposure (hereafter long and short, respec- 
tively) for the cluster OGLE-LMC0531 as an example. It is 
clear that our photometry is accurate, since bright stars lie 
very near the zero residual line. Using this diagram we can 
determine the long saturation limit, where the residuals be- 
come negative. We can also determine the faint magnitude 
limit where the short signal to noise ratio (S/N) becomes 
small and the scatter in the residuals becomes larger. A sys- 
tematic photometric offset of a few hundredths of magnitude 
where found between short and long and were interpreted 
as due to small seeing variations. They were dealt with by 
means of aperture corrections applied on the short magni- 
tudes. 

After photometric homogenization, we merge both long 
and short photometry tables. In essence, the merger process 
keeps only short stars brighter than the long saturation limit 
and only long stars fainter than the short low S/N limit. 
In the region between the quoted limits we choose the star 
with the best magnitude measurement, i.e. the star with the 
smaller photometric error. 

In Fig. [3] we show the result of the merging process. 
We also display the long saturation and short faint magni- 
tude Umits computed as previously described (solid lines). 
The CMDs display a characteristic shape, being dominated 
by the LMC fields stars, as expected, since our targets are 
intermediate to low-mass clusters. Moreover, we can distin- 
guish the old and young population of the LMC divided by 
the extended Hertzsprung-gap, and we can identify a very 
clear Red Giant Branch (RGB) followed by the Red Clump 
(RC) and the Asymptotic Giant Branch (AGB). On the blue 
side of the CMDs the Main-Sequence (MS) extend to the do- 
main of bright magnitudes, showing the presence of younger 
populations. 
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Figure 2. The residual diagrams for OGLE-LMC0531 in the B 
(lower panel), V(middle panel), and I(upper panel), where we 
show the difference between the instrumental magnitudes of the 
long and short images as a function of the long magnitude. 



Here we adopt a point spread function (PSF) fit pho- 
tometry, the standard approach when dealing with crowded 
fields. The PSF model was computed using the images from 
the standard fields in order to avoid crowding effects on 
the PSF modeling. The photometry was carried out using 
DAOPHOT code from Stetson (1994). The photometry pro- 
cess took place as following: (i) For each cluster we chose the 
deepest and cleaner field for reference, usually the long ex- 
posure image on the I filter. Automatic detection of sources 
was carried out on this reference field using DAOFIND task, 



3 ANALYSIS 

Since the clusters here studied are poorly populated objects, 
as shown in Fig. [T] it is crucial to differentiate the star clus- 
ter population from the field population. In fact, the field 
decontamination is the most challenging step in our anal- 
ysis. For instance, KMK88-38 appears on the image as an 
excess of only ~ 13 % relative to the same area in the field. 

One popular method of decontamination is the statis- 
tical com parison of CMDs such as the one developed by 
iKerber et akl (|2002l ) . In the literature other methods are 
employed with succ es s such as the st a tistic al methods of 
iGallart et al] (|2003l '): iBonatto fc Bical (|2007l ):J Piatti et al 



('2008') and the proper motion analysis made bv lRicher et al 
(,200a) 



3.1 Decontamination 



The statistical method of iKerber et all (|2002l ) is based on 
the fact that the cluster population should appear as an over 
density in the colour-magnitude space (CMS) when com- 
pared against field stars. The decontamination is made by 
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Figure 3. V,(B-V) (right panel) and V,(V-I) (left panel) CMDs from our sample. The combined long and short CMD for the entire SOI 
field is shown. The x are the stars that come from the short photometry, and the black dots are from the long long photometry. The 
upper solid line indicates the saturation limit and the lower line indicates the faint magnitudes limit. The name of the cluster associated 
to the field is indicated in each panel. 



dividing the CMD into bins of colour and magnitude and 
then comparing star counts in each bin in the cluster region 
with the corresponding bin in the field region. This method 
deals with very sharp boundaries when counting stars, which 
leads to some shot noise when working with small numbers. 
The method uses a random process to select whether each 
star is a cluster member or not. This also leads to fluctua- 
tions on the final decontamination results. 

An obvious improvement is to add an extra colour 
axis to the statistical field subtraction process. In this 
work we explore two colours simultaneously and develop 
a method of decontamination that operates on the colour- 
colour-magnitudc space (CCMS). To avoid sharp boundaries 
in the CCMS bins, each star is replaced by a 3 dimensional 
Gaussian distribution with standard deviation equal to its 
photometric error. These additional features bring all the 
information available from the photometry into the analy- 
sis, using not only the position of the star in the CCMS but 
also its associated photometric error. 



The cluster region and center are determined by eye. 
Typical values for the cluster region radii are of ~ 20" which 
corresponds to a fraction of ~ 1/60 of the SOAR/SOI FOV. 
In general the clusters are too sparse and low density to 
allow a more systematic approach. We use all the remaining 
SOI area as the field region. For both the cluster and field 
regions we take the following steps: (i) The CCMS is divided 
in a very fine grid of 0.01 in magnitude and 0.005 in both 
colours, (ii) For each bin we compute the contribution of 
each star by integrating its 3d-Gaussiaii over the bin volume, 
(iii) The resulting star counts in each bin are normalized by 
the different region areas. 

After applying step (iii) we subtract the field star counts 
from the cluster ones at each bin. We then smooth out the 
resulting decontaminated cluster counts on CCMS by re- 
binning them on a coarser grid. For visualization purposes, 
we project the coarse CCMS grid onto the usual (V-I),V and 
(B-V),V CMD planes. The coarse grid size is chosen as the 
smallest possible that is able to evidence a population that 



resembles a SSP with minimal noise and acceptable resolu- 
tion in the CMD plane. As we replaced a point process by 
smoothed 3d-Gaussians, the resulting CMDs are expressed 
in terms of star counts per bin. 

3.2 Controlled experiments 

In order to test the developed decontamination algorithm 
we performed controlled experiments using artificially gen- 
erated CMDs with parameters similar to those expected for 
the clusters in our sample. 

To generate an artificial CMD, we used the aforemen- 
tioned software developed by iKerber et al] lj2002h . This al- 
gorithm takes into account realistic photometric errors and 
the effects of CMD broadening due to unresolved binaries; 
the adopted binary fraction adopted here is 50%. 

We generated artificial CMDs for two isocrones models: 
one with log{age) = 9.1 and Z — 0.010 for which 50 stars 
were created; the other with log{age) = 8.4 and Z — 0.010 
and containing 30 stars. 15 realizations of each model were 
created. Each artificial SSP realization is inserted into the 
real images at randomly chosen positions. By doing this we 
take into account effects caused by the stochastic nature of 
the artificial SSPs and the field. 

In Figure|4]we show the resulting decontaminated V,(V- 
I) CMDs for the first 8 of the 15 controlled experiments from 
the first isocrone mentioned above {log{age) = 9.1 and Z = 
0.010, solid line). The artificial stars generated are shown as 
crosses in all panels. The mean contrast relative to the field 
is 17%. The parameters for this model are typical for the 
expected oldest clusters in our sample and should reflect the 
quality of the decontamination algorithm for intermediate 
age SSPs. We show only 8 realizations of this experiment 
for simplicity. 

Figure [5] shows the results of the controlled experi- 
ment for the other isochrone model, with a younger age 
and same metallicity. The isochrone and artiflcial stars are 
again shown in all 8 panels. The adopted age and metal- 
licity {log{age) — 8.4 and Z = 0.010) are consistent with 
the youngest clusters in our sample. To test the method at 
yet lower constrast levels, only 30 stars where generated, 
yielding a mean contrast of 11%. 

The algorithm is successfuU in recovering the enhanced 
density contrast in the CCMS regions associated to the in- 
put SSP. Both the input Main Sequence Turn-Off (MSTO) 
and the RC positions can be identifled in almost all real- 
izations of the decontamination process applied to the first 
isochrone model (Figure 2)). The extended MS locus associ- 
ated to the young SSP is also clearly visible in all but one 
of the panels of Figure (5) The dense contaminating field and 
its associated fiuctuations, however, often leave extra over- 
dense CMD locci in some cases. In particular, in a couple 
of the realizations shown, the CMD is dominated by a faint 
and red plume of stars. This is likely caused by large and 
unaccounted for photometric errors at these faint magni- 
tudes and red colours. Considering that this plume is not 
associated to any reasonable SSP at the LMC distance, it 
can be easily discarded. More troubling is the residual field 
SGB/RGB overdensity present in a couple of realizations 
shown in Figure [S] This could make the age estimate an 
ambiguous task. On the other hand, addition of an extra 
colour, as is our case, helps remove this potential ambiguity. 
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Figure 4. V, {V — I) CMDs resulting from controlled experi- 
ments using an intermediate age SSP witii 50 stars. The employed 
isochrone is shown as the solid line in all panels. Its age is 1.25 
Gyr and its metallicity is half the solar value. The artificially gen- 
erated stars are shown as crosses. The field subtracted star counts 
are coded with a grey scale, which is shown as a colour-bar on 
the right side. The bin sizes in color and magnitude are 0.2 and 
0.5 mag respectively. 
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Figure 5. V, {V — I) CMDs resulting from controlled experiments 
using a young SSP with 30 stars. The employed isochrone is shown 
in all panels as the solid line. It has an age of 250 Myrs and a 
metallicity of half the solar value. The artificially generated stars 
are shown as the crosses. The bin size and grey scale represen- 
tation of the field subtracted star counts are the same as in the 
previous figure. 

Also, in the vast majority of the realizations in both mod- 
els it is possible to identify the input SSP as the dominant 
excess of stars in the CMDs, to which isochrones may be 
visually overlaid in order to estimate their parameters. 

3.3 Parameter determination 

CMDs are widely used as a tool for parameter determi- 
nation in the liter ature (jKerber. Santiago, fc Brocatdl2007l : 
iPiatti et ahlbOOTh . although there are alternative methods 
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that can be employed in order to obtain similar results. 
These methods include integrated sp ectroscopy such as in 
lAhumada et all (|2009'); Santos ot al. (200^) and integrated 
nP 



photometry as inlHuntcr ct al., (,2003. '): Pessev et al.1 (|2008l ). 

Figures |6] to |8] show the CMDs that result from the de- 
contamination process applied to our sample clusters. They 
can be used to extract physical parameters for the clus- 
ters. The target clusters are poorly populated, having typi- 
cally from 20 to 150 stars. This limits the accuracy of clus- 
ter parameters determination. Ideally, a statistical method 
based on CMD star counts should be used to infer ages, 
metallicities, distances and redenning. This has been done 
by several of the authors quoted in this paper. However, 
in the current case, we showed in the previous subsection 
that the fluctuation in the background field is the domi- 
nant source of noise in the decontaminated CMDs. As we 
cannot single out and eliminate a prion the contaminat- 
ing field stars present in the cluster direction, the use of 
statistical methods for CMD comparison will not provide 
a significant advantage over isochrone fitting, as neither is 
capable of eliminating the effect of field contamination on 
the defined clu s ter re g ion. This is a different situation from 
iBalbinot etahl (|2009l ): lKerber. Santiago, fc Brocatj (|200'if ). 
for instance, who used statistical CMD comparisons to in- 
fer parameters and their uncertainties for rich star clusters. 
Therefore, a visual isochrone fit is adequate to estimate clus- 
ter parameters in the present c ase. Isochrones come from the 
stellar evolutionary models bylCirardi et al.1 l|2002l 'FI. 

We chose Padova isochrones although we are aware that 
systematic differences in the parameters may arise from the 
choice of one model or an other. This issue ha s been inves- 
tigated in more detail bv iKerber fc Santiagol (|2009l ). who 
compared different stellar evolution models applied to LMC 
clusters. They show t hat, in relation to Pado va isochrones, 
the ones from BaSTI (|Pietrinferni et al] 2004) tend to over - 
estimate the age by ~ Q.ZGvrs. PEL (|Castellani et al.ll2003l ') 
isochrones may lead to a discrepancy in the distance of 
~ 2.8A:pc at SOfcpc. The effects of convective overshooting 
(CO) are analysed in the same work. Obvious differences 
in age arise between CO models and non-CO models, in 
the sense that the former yield an older age when com- 
pared to the latter (loq(ac] e),co)-loq{age)ino n -c(n - 0.1). 
Apart from that, however. iKerber fc Santiagol (l2009l ) see no 
other significant trend in the parameters. The authors also 
compared their met allicity values with th o se obtained from 
spect roscopic data (jOlszewski et al.lll99ll : iGrocholski et all 
I2OO6I ) and found that the Padova models give the best agree- 
ment between spectroscopic and photometric data. We thus 
conclude that, for the purpose of this paper, the influence 
of the model adopted is negligible and the differences from 
model to model are much smaller than the uncertainty of 
our parameter determination. 

For each cluster we derive the following parameters: age, 
metallicity, distance and colour excess. The parameters must 
obey some constraints in metallicity and distance, which 
are typical of LMC clusters. The allowed ranges adopted 
here are: 18.25 ^ (m - A/)o 18.75 and \Zq Z ^ Zq 
l|Kerber. Santiago, fc Br ocato 2 0071). We take th e reddening 
tables from Schleeel. Finkbeiner. fc David l|l998l ) as a initial 
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guess of E{B — V). In some cases, the lack of clear CMD 
structures, such as the MSTO, RGB and RC prevent strong 
constraints on the parameters, specially the distance modu- 
lus. In such cases we keep it flxed at the mean LMC value. 

The results of our parameter determination are listed in 
Table [21 We will now comment the results for each cluster 
set. In Fig. [6]to[8]we show, along with the subtracted CMDs, 
three isochrones overploted [solid lines). All isocrones have 
the same metallicity and only differ by age, whose values are 
given on the figures. On all figures the V,(B-V) CMD is on 
the left side and the V,(V-I) is on the rigth side. The cluster 
names and other isochrone fit parameters are also shown. All 
CMDs show star count in bins of colour-magnitude instead 
of stars itself, the grey scale is adjusted to minimize the noise 
and increase the contrast for better visualization. 

3.3.1 OGLE-LMC0531 and OGLE-LMC0523 

These two clusters are on the same SOAR/SOI field. Only 
the first one is an age-gap candidate. As seen in Fig. [6] the 
two clusters have markedly distinct ages. In fact, OGLE- 
LMC0531 is older, with an age between ~ 1 — 2Gyrs. No 
sign of an upper MS is seen in either the V ,{B — V) or 
V,{V - /)CMDs. In the former, its MSTO and lower MS 
are more clearly defined, along with a trace of its RC. As ex- 
plained earlier, this allowed better constraints on its physical 
parameters, yielding a E{B - V) ~ 0.09, (m - M)o = 18.35 
and a. Z = 0.013. The exti nction value is very similar to 
the initial guess taken from ISchlegel. Finkbeiner. fc DavisI 
(1998). 

OGLE-LMC0523 has an age of - 2mMyrs. In reality, 
this is an upper limit, since no clear turn-off is seen up to 
the saturation limits. Since no other clear CMD structure 
besides the MS is detected, it was harder to constrain the 
cluster distance. We adopted the mean LMC distance of 
(m - M)o = 18.50 and found Z = 0.010 and E{B -V) ^ 
0.06. 

3.3.2 KMK88-38, KMK88-39 and OGLE-LMG0214 

KMK88-38 is our original target in this field, as an age- 
gap candidate. The other clusters present on the field are 
KMK88-39 and OGLE-LMC0214. This latter was imaged 
only in the V and I bands, due to variations in telescope 
pointing. Even though we do not show its CMDs, it has 
been submitted to the same analysis as the other clusters, 
the results of which are shown in Table [2] 

The field subtracted CMDs for KMK88-38 and KMK88- 
39 are shown in Fig. [7] The panels and symbol conventions 
are the same as in Fig. [6] KMK88-38 is relatively old, with 
an age between 1 ~ 2Gyrs, similar to OGLE-LMC0531. The 
presence of strong CMD features, such as a MSTO, and an 
RGB allowed us to explore a wider range of parameters with 
the effect of better constraining their determination. The 
derived parameters for KMK88-38 are: (m — M)q = 18.35, 
Z = 0.005 and E{B ~ V) = 0.02. 

On the other hand, KMK88-39 is very young, with an 
age of ~ 250Myrs when using a metallicity oi Z = 0.009. 
The age derived for this cluster is again an upper limit, since 
no clear MSTO is seen. The distance adopted was again the 
mean LMC distance, and we found a redenning of E(B — 
V) = 0.01. 
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Figure 6. Panel a: V, (B - V) CMD of OGLE-LMC0531; panel 
6: V, {V — I) CMD for the same cluster as in previous panel. The 
same model, whose isochrone parameters are indicated, was used 
on both CMDs. Panel c: V, (B - V) CMD of OGLE-LMC0523; 
panel it d: V, {V — I) CMD for the same cluster as in panel c. The 
same model, whose isochrone parameters are indicated, were used 
on both CMDs. In all panels, the field subtracted star counts are 
coded with a grey scale, which is shown as a colour-bar on the 
right side. The adopted bin size is 0.55 in magnitude and 0.30 in 
colour. The dashed line in the upper panels marks the position of 
the MSTO of a cluster with 3Gyrs and half solar metallicity. 



3.3.3 NGC 1878 

NGC 1878 was accidentally observed in place of the age gap 
candidate cluster BSDL917. In contrast to the other fields, 
only one cluster was imaged this time. Even though is an 
NGC object, NGC 1878 had never been studied before. 

Since NGC 1878 is a relatively rich cluster (> 100 stars) 
it allowed us to explore our decontamination algorithm in a 
simpler case, where the density of cluster stars in the CMD 
at a given colour-magnitude bin is significantly greater than 
the noise generated by the fluctuation of field stars in this 
same bin. In Fig. [8] we show its subtracted CMD, adopting 
the same panels and symbols as in the two previous figures. 
We clearly see much more well defined CMD loci than in the 
other cases, indicating that the decontamination process is 
effective. 

The V,(B-V) CMD of NGC 1878 looks slightly bluer 
than expected from the V,(V-I) CMD. This means that 
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Figure 7. Same as in Fig.[6]but for the clusters KMK88-38 (pan- 
els a and b) and KMK88-39 (panels c and d). The adopted bin 
size is 0.40 in magnitude and 0.25 in colour 



Name 


log{Age) 


Z 


E{B - V) 


(m - M)o 


OGLE-LMC0214 


8.4 i 


0.013 


0.10 


18.50* 


OGLE-LMC0523 


8.3 4. 


0.010 


0.06 


18.50* 


OGLE-LMC0531 


9.1 ±0.1 


0.013 


0.09 


18.35 


KMK88-38 


9.1 ±0.1 


0.005 


0.02 


18.35 


KMK88-39 


8.4 i 


0.009 


0.01 


18.50* 


NGC1878 


8.3 ±0.2 


0.009 


0.12 


18.50* 



Table 2. The result of our CMD analysis. We list the object 
name, the derived log(age), metalicity, E{B — V) and the distance 
modulus. The arrows indicate upper limits. The asterisk indicates 
where the distance modulus was fixed at the given value. 



a unique isochrone solution to both CMDs was harder to 
achieve. As a compromise, the best visual fit solution is 
slightly red (blue) when compared to the main MS locus in 
the former (latter) CMD. This effect can result from a com- 
bination of crowding and variable PSF, which may lead to 
increased light contamination on PSF fitting to the B band, 
therefore making the star seem bluer. The best isochrone fit 
gives a young age, ^ 200Myrs, a typical LMC metaUicity 
oi Z = 0.009, E{B - V) = 0.12. We adopted a distance 
modulus of (m - M)o = 18.50. 
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Figure 8. Same as in Fig.[6]but for the cluster NGC 1878. The 
adopted bin size is 0.50 in magnitude and 0.25 in colour 



4 CONCLUSIONS 

We have obtained SOAR/SOI images of 3 fields containing 
6 LMC clusters, two of which were previously identified as 
candidates to fill the LMC age-gap. Photometry was car- 
ried out in three filters (B,V and I) reaching well below the 
MSTO of an old stellar population. 

All but one of these clusters are poor and sparse, con- 
sistent with intermediate to low initial masses, requiring a 
careful and objective process to eliminate field stars from the 
cluster CMD. We developed such method that simultane- 
ously uses the information in the available three-dimensional 
space of magnitudes and colours. With the resulting decon- 
taminated CMDs we carried out visual isochrone fits based 
on Padova evolutionary sequences. 

Our main result is that none of the clusters here stud- 
ied is older than 2 Gyrs, therefore not filling the LMC 
age-gap. The two main candidates, OGLE-LMC0531 and 
KMK88-38, are in fact the oldest in the sample, display- 
ing a clear MSTO and a branch of evolved stars. Their 
MSTO magnitudes are about 2 mag brighter than that ex- 
pected for a 3 Gyrs old SSP at the LMC distance, as evi- 
denced from Figures [6] and [7] Interestingly, they are sided 
on the images by much younger clusters, for which only up- 
per limits in age (< 200Mj/rs) could be derived from the 
CMDs. This large age range is evidence either of a com- 
plex cluster formation history, or very efficient orbital mix- 
ing within the LMC. Another noticeable result is that the 
ages here inferred are similar to the last two close encoun- 
ters between the LMC and SMC according to rec ent N-body 
simulations (jCardiner. Sawa. fc Fuiimotd Il994l ) , adding a 
few more examples of clusters that likely have formed 
as the result of SMC/LMC close- by passages. We also 
note that the two older clusters are systematically on the 
foreground relative to current models for the LMC d isk 
l|Kerber, Santiago, fc Brocatcll2007l : iNikolaev et"al]|2004) . 

The clusters in our sample are of lower initial mass 
(M ~ 10* Mq) than most LMC clusters studied so far. The 
fact that they do not fill the age gap may indicate that yet 
lower mass clusters, possibly cluster remnants, must be sam- 
pled and studied in order to bridge the apparent inconsis- 



tency between reconstructed cluster and field star formation 
histories in the LMC. 

Acknowlegments.This work was supported by Con- 
selho Nacional de Desenvolvimento Cienti'fico e Tecnologico 
(CNPq) in Brazil and Coordenagao de Aperfeigoamento de 
Pessoal de Nfvel Superior. 



REFERENCES 

Ahumada A. V., Talavera M. L., Claria J. J., Santos 
J. F. C, Bica E., Parisi M. C, Torres M. C, 2009, lAUS, 
256, 293 

Balbinot E., Santiago B. X., Bica E., Bonatto C, 2009, 
MNRAS, 396, 1596 

Baume C, Noel N. E. D., Costa E., Carraro G., Mendez 
R. A., Pedreros M. H., 2008, MNRAS, 390, 1683 

Bica E. L. D., Schmitt H. R., Dutra C. M., Oliveira H. L., 
1999, AJ, 117, 238 

Bonatto C, Bica E., 2007, MNRAS, 377, 1301 

Castellani V., Degl'Innocenti S., Marconi M., Prada Mo- 
roni P. G., Sestito P., 2003, A&A, 404, 645 

de Grijs R., Anders P., 2006, MNRAS, 366, 295 

Gallart C, Zoccali M., Bertelli G., Chiosi C, Demarque 
P., Girardi L., Nasi E., Woo J.-H., Yi S., 2003, AJ, 125, 
742 

Gardiner L. T., Sawa T., Fujimoto M., 1994, MNRAS, 266, 
567 

Girardi L., Bertelli G., Bressan A., Chiosi C, Groenewegen 
M. A. T., Marigo P., Salasnich B., Weiss A., 2002, A&A, 
391, 195 

Grocholski A. J., Cole A. A., Sarajedini A., Geisler D., 

Smith V. V. 2006, AJ, 132, 1630 
Holtzman J. A., et al., 1999, AJ, 118, 2262 
Hunter D. A., Elmegreen B. G., Dupuy T. J., Mortonson 

M., 2003, AJ, 126, 1836 
Javiel S. C, Santiago B. X., Kerber L. O., 2005, A&A, 431, 

73 

Jensen J., Mould J., Reid N., 1988, ApJS, 67, 77 
Kerber L. O., Santiago B. X., 2009, lAUS, 256, 391 
Kerber L. O., Santiago B. X., Brocato E., 2007, A&A, 462, 
139 

Kerber L. O., Santiago B. X., Castro R., Valls-Gabaud D., 

2002, A&A, 390, 121 
Kontizas E., Metaxa M., Kontizas M., 1988, AJ, 96, 1625 
Nikolaev S., Drake A. J., Keller S. C, Cook K. H., Dalai 

N., Griest K., Welch D. L., Kanbur S. M., 2004, ApJ, 601, 

260 

Noel N. E. D., Gallart C, Costa E., Mendez R. A., 2007, 
AJ, 133, 2037 

Olszewski E. W., Schommer R. A., Suntzeff B., Harris H. 

1991, AJ, 101, 515 
Parmentier G., de Grijs R., 2008, MNRAS, 383, 1103 
Pessev P., Goudfrooij P., Puzia T., Chandar R., 2008, AAS, 

211, #162.27 

Piatti A. E., Geisler D., Sarajedini A., Gallart C, 2009, 

A&A, 501, 585 
Piatti A. E., Geisler D., Sarajedini A., Gallart C, Wischn- 

jewsky M., 2008, MNRAS, 389, 429 

Piatti A. E., Sarajedini A., Geisler D., Gallart C, Wischn- 
jewsky M., 2007, MNRAS, 382, 1203 



LMC clusters at intermediate masses 9 

Pietrinferni A., Cassisi S., Salaris M., Castelli F., 2004, 
ApJ, 612, 168 

Pietrzynski G., Udalski A., Kubiak M., Szymanski M., 
Wozniak P., Zebrun K., 1999, AcA, 49, 521 
Rich R. M., Shara M. M., Zurck D., 2001, AJ, 122, 842 
Richer H. B., et al., 2008, AJ, 135, 2141 
Santiago B. X., 2009, lAUS, 256, 69 

Santos J. F. C, Jr., Claria J. J., Ahumada A. V., Bica E., 
Piatti A. E., Parisi M. C, 2006, A&A, 448, 1023 
Sarajcdini A., 1998, AJ, 116, 738 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Sharpcc B., Stark M., Pritzl B., Smith H., Silbermann N., 

Wilhelm R., Walker A., 2002, AJ, 123, 3216 
Stetson P. B., 1994, PASP, 106, 250 



